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Abstract 

We study the so-called "sunset diagram", which is one of two-loop self-energy diagrams, for 
scalar field theories at finite temperature. 

For this purpose, we first find the complete expression of the bubble diagram, the one-loop 
subdiagram of the sunset diagram, for arbitrary momentum. The temperature-dependent discon- 
tinuous part as well as the well-known temperature independent part is obtained analytically. The 
continuous part is reduced to a one-dimensional integral which one can easily evaluate numerically. 

We calculate the temperature independent part and dependent part of the sunset diagram sep- 
arately. For the former, we obtain the discontinuous part first and the finite continuous part next 
using a twice-subtracted dispersion relation. For the latter, we express it as a one-dimensional 
integral in terms of the bubble diagram. 

We also study the structure of the discontinuous part of the sunset diagram. Physical processes, 
which are responsible for it, are identified. Processes due to the scattering with particles in the 
heat bath exist only at finite temperature and generate discontinuity for arbitrary momentum, 
which is a remarkable feature of the two-loop diagrams at finite temperature. 

As an application of our result, we study the effect of the diagram on the spectral function of 
the sigma meson at finite temperature in the linear sigma model, which was obtained at one-loop 
order previously. At high temperature where the decay a — > irir is forbidden, sigma acquires 
a finite width of the order of lOMeV while within the one-loop calculation its width vanishes. 
At low temperature, the spectrum does not deviate much from that at one-loop order. Possible 
consequences with including other two-loop diagrams are discussed. 
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I. INTRODUCTION 



The so-called sunset diagram, which is depicted in Fig.l, is one of two- loop self-energy 
diagrams. The sunset diagram appears in theories with 4-point vertices such as 0(4) linear 
sigma model, 4 theory, and so on. 

At finite temperature, two-loop self-energy diagrams have remarkable features that are 
not seen at zero temperature. Consider the discontinuity of the self-energy. At zero temper- 
ature, the discontinuity of one-loop diagrams is due to two particle intermediate states and 
that of two-loop diagrams comes from three particle states in addition to two particle ones. 
In general, in higher loop diagrams new processes appear. However, they contribute to the 
discontinuity only at higher energies. On the other hand, at finite temperature, even though 
the number of particles participating in the process at a given order of loops is the same as 
that at zero temperature, new processes appear even at low energy. This is because at finite 
temperature some of particles participating in the process can be particles in the heat bath. 
Furthermore, as will be shown in this paper, there exist processes which are possible at 
arbitrary energy. Accordingly, the discontinuity of two-loop self-energies is non-vanishing in 
all the energy region. Therefore, in some cases at finite temperature, to extend calculations 
to two-loop order has a meaning more than just making more precision. 

One of such cases is the spectral function for the sigma meson in the linear sigma model. 
The spectral function of the sigma meson at finite temperature was studied by Chiku and 
Hatsuda at one- loop level £Q. The sigma meson at zero temperature has a large width due 
to the strong coupling with two pions. However, they found that at finite temperature the 
spectral function near a — ► mi threshold is enhanced as a typical signal of chiral phase 
transition. This is because, as the temperature increases, the mass of a decreases while that 



FIG. 1: "Sunset" diagram. We label the external particle by "<E>" and internal particles with mass 
rrii by <fo. 
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of 7r increases due to partial restoration of chiral symmetry and, accordingly, the phase space 
available for the a — > decay is squeezed to zero. At finite temperature, however, there 
exist processes, collision with or absorption of a thermal particle in the heat bath, which 
contributes to the discontinuity at arbitrary energy for a, as discussed above. If we include 
them, the structure of the spectral function might be significantly modified. Their effects on 
the spectral function cannot be taken into account till we extend the calculation including 
two-loop self-energies. 

In this paper, we take a step forward in the calculation of two-loop self-energy diagrams, 
i.e. we evaluate the thermal sunset diagram for scalar field theories. After giving a brief 
review of the real time formalism [2TTHj in the second section, which we use for the calculation 
of diagrams, we first examine the bubble diagram, a one-loop diagram which appears as a 
subdiagram of the sunset diagram, at finite temperature in the third section. Then, we 
discuss the structure of the discontinuous part of the thermal sunset diagram in the fourth 
section and explain how we calculate the sunset diagram in the fifth section. Using our 
result for this diagram, we study the effects of the diagram on the spectral function of a at 
finite temperature in 0(4) linear sigma model in the sixth section. Finally we summarize 
the paper in the seventh section. 

II. BRIEF REVIEW OF THE REAL TIME FORMALISM 

For the calculations of thermal Feynman diagrams we adopt the real time formalism 
throughout this paper. We briefly review the formalism in this section. 

In the real time formalism propagators are given by 2 x 2 matrices. The 4-components 
of the free propagator of a scalar particle with mass m are given by ^3] 




(II.l) 



(II.2) 




(II.3) 




—i 




where n(k Q ) is the Bose-Einstein distribution function at temperature T = 1/13: 



Off-diagonal elements depend on a free parameter a. In this paper we make the symmetrical 
choice a = (3/2, leading to 

tA F 2 (k; m) = iA^k; m) = e^ /2 n(k )2^5(k 2 - m 2 ). (II.6) 

We denote the matrix of the full propagator by 

A( fc ; m )=f An( * ; " l) ^'"'V (II.7) 
yA 2 i(k;m) A 22 (k;m) J 

Each component is the solution of the Schwinger- Dyson equation: 

A ab (k; m) = A F b (k; m) + A F c (k; m)U cd (k)A db (k; m), (11.8) 

where U cd (k) is the self-energy 
The following matrix 



U(k) = I >/!+»(*») ^ I , (II.9) 
\/n(k Q ) ^l + n(k ) 

'diagonalizes' the free propagator: 

Af J k;m) A F (k;m)\ , ( A F (k;m) . 

U(k) ( 111 ' J 12V ' ; I U(k) = 01 ' , (11.10) 

A£(*;m) A F 2 (k;m) J \ -A F (k;m) 

where A F (k; m) is the Feynman propagator at T = 

A o(k;m) = - L— . (11.11) 

k l — m l + it] 

It is known that U(k) also diagonalizes the full propagator ^3] as well as U a b{k) 

. n(k) . 

U ab (k) = (U(k) )ae (U(k) ) db . (11.12) 

o -mrj cd 

This matrix equation gives relations between matrix elements: 

Refl(fc) = Ren n (fc), (11.13) 
Imfl(A;) = tanh(/3|A;o|/2)Imn 1 i(A;), (11.14) 
lmR{k) = isinh(/?|A; |/2)n 12 (A;). (11.15) 
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Let us next find the expression of the spectral function. An has the following spectral 
representation: 

A u (k;m) = J dcu 2 p(uj,k)A( 1 (k ,uj), (11.16) 

where we denote the argument of A n by /c and to = \/PTff? instead of k and m. One 
can prove that, if An is given in Eq. (|II.16j) . the Green function 

A(t,x) = j^yf d 4 ke ik -*A n (k;m) (11.17) 

obeys the Kubo-Martin-Schwinger (KMS) condition ^3]. From Eq. ()11.16|) we obtain 

p(k) = — tanh(/5|A;o|/2)ImA 11 (A;; m). (11.18) 

7T 

Our next task is thus to find the expression of ImAn(A;; m). Denoting the diagonalized full 
propagator by 

I A(jfe;m) , 
U{k)- 1 A{k- 1 m)U{k)- 1 = \ V ; , (11.19) 

y -A(k;m) 

we obtain 

A(k~ m) 
A = U(k) | | U(k) 

-A(k;m) 



[1 + n(fc )]A(fc; m) - n(k )A(k; m)* y/n(ko)[l + n(k )}[A(k; m) - A(k; m)*\ 

^/n(k )[l + n(k )}[A(k; m) - A(k; m)*\ -[1 + n(k )]A(k; m)* + n(k )A(k; m) 

(11.20) 

Taking the (1, l)-components in the both sides of the above equation yields 

ImAii(Jfe; m) = coth(/3|jfe |/2)IniA(ife; m). (11.21) 



The expression for A(k;m) can be found in the following manner. In Eq. (jII.8|) . iteratively 
using Eq. (jII.8j) itself in the right hand side and diagonalizing, we obtain 

A^(k;m) 

-A^(k;m)* 



U(k)- l A(k;m)U(ky l 



f A^(k;m) \ / U(k) \ / A^(k;m) ( 

-A^(k;m)* ) \ -5(fc)* / \ -£$(k;m) 



(11.22) 
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The (1,1) component of this equation gives 

A(k;m) = (f/(A;)- 1 A(A;;m)f/(A;)- 1 ) 11 

= Aq (k; m) + Aq (k; m)fl(k)A% (k;m)-\ 

= k 2 - m 2 - Wk) ' ^ IL23 ^ 
Therefore, from Eqs. (jII.18|) . f)II.21[) and (jll.23|) . we obtain the desired expression of the 
spectral function: 

1 Imll(fc) 
P ^ ' ~ ~n [k 2 - m 2 - Rell(fc)] 2 + [ImIl(A;)] 2 ' ( 2 ' 

Thus the spectral function can be written only with one component of the self-energies while 

all the components enter A n (A;;m) (see Eq. (jll.8p ). 

III. BUBBLE DIAGRAM 

In this section we examine the "bubble" diagram, which is a one-loop diagram shown in 
Fig. 2. The bubble diagram appears as a subdiagram of the sunset diagram and therefore 
we need the expression of the former in the calculation of the latter. 

4>i 




FIG. 2: "Bubble" diagram. 

Among the 4-components of the bubble diagram we need only the (1,1) component since 
the sunset diagram does not have any internal vertices. The (1,1) component of the bubble 
diagram for scalar particles is given by 

/d 4 k 
-^—rgiA 11 (p + A;;mi)iAf 1 (A;;m2). (III. 25) 

Eq. (|111.25|) can be expressed as the sum of terms with different numbers of Bose-Einstein 
factors, n: 

;mi,m 2 ) + (F^ip; m 1 ,m 2 ) + (l <-> 2)) + F^(p; mi,m 2 ). (III. 26) 



Each term in the right hand side is respectively given by 



[p ;m ll m 2 ) 
;mi,m 2 ) 
;m 1 ,m 2 J 



(27r) 4 (p + A;) 2 — m\ + in k 2 — m 2 , + irj' 



d 4 k 



-n{k )2nS{k 2 - m\) 



[111.21) 
;ill.28) 



(27r) 4 (p + A;) 2 -mj + ir] 

-^rjnipo + k )2n5((p + kf - mj)n(k )2n5(k 2 - m 2 ). (111.29) 
(2tt) 4 



Let us carry out the integration in the above equations. The T-independent part 1^ is 
given in textbooks |16j . In <i-dimension it is given by 

,2 ^2 



7 (2) (p 2 ;mi,m 2 ) 



16tt 2 



1 , m; , m 9 
-- + x+ln^- - x_ln^f - 2-1 



K 



where k is the renormalization point and 4 and x± are 
1 1 



e e 



— 7 + ln47r (e = (4 — d)/2, 7 : Euler constant) 



1 2 2 

1 mf — 



2p 2 



(111.30) 



(111.31) 



(111.32) 



/ is given by 



In 



(x + -VC)(x-+VC) 



+ i2ir 



-2y/D 



arctan-^t — arctan4= 

V D v D 



Cln 



where C and -D are 



(v / C-x + )(v / C+a:_) 
(VC-x^)(VC+x+) 



for p 2 > (mi + m 2 ) 2 , 

for (mi — m 2 ) 2 < p 2 < (mi + m 2 ) 2 , 

for p 2 < (mi — m 2 ) 2 , 



C =4 



(mi + m 2 ) 
p2 



(rrii — m 2Y 
p2 



D = -C. 



(111.33) 



(111.34) 



We see from Eas. fllll.27jl and (1111.331) that the discont inuous part, Imi/( 2 \ is non- vanishing 
for p 2 > (mi + m 2 ) 2 . 

The divergent part in Eq. (jIII.27j) is renormalized by applying a counter term to the 
Lagrangian. The finite part is determined in such a way that the resultant self-energy 
satisfies a proper normalization condition. 

Let us next turn to the T-dependent part, and F^ 3 \ We first discuss their discon- 
tinuous parts. As shown in Appendix 1X1 the discontinuous part of F' 2 '(p; mi, m 2 ) is given 
analytically as follows: 
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1. For p 2 > (mi + nT-2) 2 or < p 2 < (m x — m 2 ) 2 , 



lmiF^ 2 \p; mi, m 2 ) 



1 1 



In 



1 - e~^+ 



16tt|p| y0 1-e-^- 

2. For (mi — m 2 ) 2 < p 2 < (mi + m 2 ) 2 , 

ImiF^(p; mi, m 2 ) = 0. 

3. For p 2 < 0, 



(111.35) 



(IIL36) 



ImiF^ (p; m l5 m 2 ) 



1 -1 
167r|p|7" 



In |(1 -e-^+Xl -e _ ^-)| . (111.37) 



Here, 



2 2\ 2 

m 2 — mf x 



bo|± 



(m 2 + mi) 2 (m 1 -m 2 ) 



p* 



p z 



(IIL38) 



i 7 ^ 3 ) has only discontinuous part. Its calculation can be done in a way similar to that for 
lmiF^(p; m 1: m 2 ). We show the final results: 



1. For p 2 > (mi + m 2 ) 2 , 

F (3) (p;mi,m 2 ) 



1 1 



-In 



I _ e -/3^+ e /3(|po|-w-) _ I 



8ir\p\(3ePM - 1 

2. For (mi — m 2 ) 2 < p 2 < (m x + m 2 ) 2 , 

F^(p;mi,m 2 ) = 

3. For < p 2 < (mi ~ ^) 2 , 

_p(3) ^. m ^ m ^ 



1 _ e -/3o;_ e /3(|p |-w+) _ 1 



(111.39) 



(111.40) 



1 


8tt 


p|/5 


X 


In 



1 - e~P"+ 



1 - e-^- 



3 /3|po| ln 



1 _ e -P{\Po\+u+) 



1 _ e -P(\ P0 \+u>-) 



(111.41) 



4. For p 2 < 0, 



F( 3 )(p;mi,m 2 ) = ^-jL_ j^-^ [-ln|l - e~^\ + e -^ln|l - e -^~^\] 
+ _J — _ [_l n ji _ e -^-| + e^llnll - e - /3( "- +|pol) |] | . 



,P\po\ 



(111.42) 



9 



In the above equations uj± are given by Ea. ()111.38jl with mi and m<i replaced by max{mi, 1712} 
and min{mi, 777,2}, respectively, since F^(p; mi, 777-2) is symmetric with respect to mi and 
777 2 by definition. 

It should be noted that the T-dependent part has two cuts in the complex p 2 plane: 
one starts from p 2 = (mi + 77i 2 ) 2 to the right along the real axis and the other from p 2 = 
(mi — 7772) 2 to the left. 

Let us next discuss the continuous part of F^ 2 \ which is 

RetF^(p;mi,m 2 ) = - [ ^-V—-^ ^n(k )27r5(k 2 - m\), (111.43) 

J (27t) 4 (p + ky — m{ 

where V stands for the prescription of Cauchy's principal value. After integration over fco 

and angle, we can express Eq. (jIII.43|) as a one-dimensional integral: 

ReiF- 2 -(p; mi,m 2 ) 

1 f°° 

= - — - / dun(uj) 

16tt 2 |p| J m2 

(p 2 + 2p uj — 2|p|-\/a; 2 — m\ + m\ — m\)(p 2 — 2p uj — 2|p|a/cg> 2 — m\ + m 2 . — m\) 



x Pin 



(p 2 + 2pou] + 2|p| \J uj 2 — m\ + 777 2 — m 2 x )(p 2 — 2p uj + 2|p| \Ju 2 — m\ + m\ — m\ 

(111.44) 



whose integration will be carried out numerically. 



IV. STRUCTURE OF THE DISCONTINUOUS PART OF THE SUNSET DIA- 
GRAM 

Before proceeding to the calculation of the thermal sunset diagram, we study the structure 
of the discontinuous part of the diagram. In ref.^H]; Weldon analyzed the discontinuous part 
of the bubble diagram. We will generalize it for the sunset diagram. For this purpose it is 
convenient to use Eq. (jII.15j) |17j . namely 

ImiX sun (A;; m x ,m 2 , m 3 ) = i sinh(/5|/c |/2)7Xsun(^; m x ,m 2 , m 3 )i 2 . (IV.45) 

Here -^"sun ( ^ ] 

mi, m 2 , 7773) is the (1,1) component of the diagonalized self-energy matrix for 
the sunset diagram. We note that the self-energy which actually enters spectral functions is 
diagonalized one (see Eq. (|II.24jl ). I sun (k; mi, m 2 , 7773)12 is the (1, 2) component of the sunset 
diagram, which is given by the following integral: 

f d A p 

il sun (k;mi,m 2 ,m 3 )i 2 = J j—^iA( 2 (p; mi)il huh (k - p; m 2 , m 3 )i 2 . (IV.46) 
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For iAf 2 (p;mi) it is convenient to use another form: 

iAf 2 (p;mi) = e^ 2 f(p )e(p )2ir5(p 2 - m\) 



/fa) = ^n- ( IV - 47 ) 

i2b u b(p; 777,1, 7712)12 denotes the (1,2) component of the bubble diagram: 

il huh (k - p;m 2 ,m 3 ) 12 = i J j^^iA( 2 (q; m 2 )iA( 2 (k - p - q; 777,3). (IV.48) 
One can reduce this equation to 

i2bub(p;777i,m 2 )i2 
= ie P{ko ~ Po)/2 f{k -p ) 

r d 4 o 1 

x / m \iaj? tp { 1 + ra 2 + ™3 <J(*o -p -E 2 - E 3 ) - 5(k -p -E 2 + E 3 )] 

J {Z7C) Z 4:E 2 E 3 

-(77 2 - 77 3 ) [5(&q - Po - E 2 + E 3 ) - 5(k - Po + E 2 - E 3 )}} , (IV.49) 

where E 2 = \/\q,\ 2 + ml and E 3 = yjk — p — q| 2 + m\. We also define E 1 = a/IpI 2 + 777f 
for later use. 77j is the Bose-Einstein factor defined by 77j = n{Ei). Using Eq. ()lV.45|) and 
Ea. (llV.4(ij) in which Ea. (llV.49l) is substituted yields 

ImzXgun^; 777,1, 777,2,7773) 

d 3 p f d 3 q 1 



(2tt) 3 J (2tt) 3 8E ± E 2 E 3 

X {((1 + m)(l + 77 2 )(1 + 77 3 ) - 77 1 77 2 77 3 )5(A:o - E x - E 2 - E 3 ) 
+ (77l77 2 77 3 - (1 + m)(l + 77 2 )(1 + 77 3 ))5(A:o + E 1 + E 2 + £3) 
+ [(77 2 77 3 (1 + Til) - 77i(l + 77 2 )(1 + Tl 3 ))5(k - E X + E 2 + £3) 

+ (771(1 + 77 2 )(1 + 773) - 77 2 77 3 (1 + 77i))5(A;o + E\ — E 2 — E 3 ) + (permutations)]} . 

(IV.50) 

In this equation "permutations" stands for the terms obtained by permuting the particle 
labels of the third and the forth terms. 

Let us now consider the physical content of Eq. (jIV.50|) . The first term in Eq. (jIV.50|) 
may be interpreted as the probability for the decay $ — > <pi<f) 2 4> 3 with the statistical weight 
(1 + 77^(1 + 77 2 )(1 + 773) for stimulated emission minus the probability for the creation 
0i<^203 —>■ $ with the weight n\n 2 n 3 for absorption. The second term is the anti-particle 
counter part of the first term. The third term represents the probability for $0 2 03 — > 0i with 
the weight 772773(1 + 77!) minus that for (pi — > §4> 2 (j) 3 with the weight n 1 (l + 772)(l + 77 3 ). Here 
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fa stands for the anti-particle of fa. The forth term is the anti-particle counter part of the 
third term. It represents the probability for $0i — > 0203 with the weight ni(l + n 2 )(l + n 3 ) 
minus that for 0203 — * $<f>i with the weight 712713(1 +n>i). All processes are shown in FigJHl 
We next find the region of k 2 where the physical processes contained in Eq. pV.50jl are 
possible, which is equivalent to looking for the condition under which the integral over q in 
Eq. (|IV.50|) survives. For the first and the second terms in Eq. flIV.50|) to be non- vanishing, k 2 
must satisfies the condition k 2 > (777.1 + M 23 ) 2 , where M23 is the invariant mass of 02 and 03. 
Therefore, the processes in Fig|3](a), (b), (c) and (d) are possible for k 2 > (m 1 + m 2 + 77i 3 ) 2 
since M23 > 7712 + 7773. The third and the forth terms survive when k 2 < (mi — M23) 2 . 
The processes in FigOl (e), (f), (g) and (h), therefore, take place at arbitrary k 2 . This is 
reasonable since the processes in FigEl (g) and (h) are scattering ones and since those in 
Figim (e) and (f ) can be also regarded as scattering by interpreting the incoming $ in (e) as 
an outgoing anti-$ and by doing the outgoing $ in (f) as an incoming anti-<3>. Accordingly, 
the discontinuous part of the sunset diagram is non- vanishing for arbitrary k 2 , which is a 
remarkable feature of the thermal self-energy at and beyond two-loop order. 



V. CALCULATION OF THE THERMAL SUNSET DIAGRAM 

In this section, we explain how to calculate the thermal sunset diagram by reducing it to 
an expression written in terms of the bubble diagram previously obtained. 
The (1,1) component of the sunset diagram shown in Figgis given by 

— -piAf^pjTOi) J ^y I zAf 1 (g;777 2 ) 7Af 1 (A;-p-g;m3). (V.51) 

We decompose Eq. (|V.51J) into terms without and with Bose-Einstein factors, which we 
denote by X^(k 2 ; m%, m 2 , m 3 )n and T^ n {k; mi, m 2 , m 3 )n, respectively: 

l sun (k; mi, m 2 , m 3 ) u = I^^(k 2 ; m a , m 2 , m 3 )n + I^ n {k; m u m 2 , m 3 ) u . (V.52) 

They are given by 

:QA; 2 ; mi, m 2 , 7773)11 = /t^^ ^— /( 2 '((^ - p) 2 ; ^2, (V.53) 

J (2vr) 4 p z — m| + vr] 



l t sm (k;m 1 ,m 2 ,m 3 ) n 



d 4 p 



x 



(27r) 4 p 2 — m 2 + irj 

[(F {2) (k - p; m 2 , m 3 ) + (2 «-> 3)) + F {3 \k - p; m 2 , m 3 )] 
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/d 4 p 
j^-^n(p )27r5(p 2 - m\) [l {2) ((k - p) 2 ;m 2 ,m 3 ) 

+ (F ( - 2 \k-p;m 2 ,m 3 ) + (2 <-> 3)) + T (3) (A; - p; m 2 , m 3 )] . 

(V.54) 

Here we have expressed the second integral in Eq. (|V.51|) in terms of I^ 2 \ and F^ 3 \ The 
T-independent part, Eq. (jV.53|) . has a subdivergence coming from the nested bubble diagram 
and a two-loop overall divergence. On the other hand, T-dependent part, Eq. ()V.54j) . has 
only a subdivergence, which can be removed by carrying out renormalization at one-loop 
level. Hereafter, 1^ expresses that with the divergence removed. We calculate Eqs. flV.53|) 
and (jV.54|) separately. 



A. T-independent part 

The T-independent part of sunset type diagrams has been calculated by several authors 
|19ll2*2*] so far. In this paper we calculate its finite part using a dispersion relation. 

We first find the discontinuous part of 2^(/c 2 ; mi, m 2 , m 3 )n and then compute the finite 
continuous part using the obtained discontinuous part, via the twice-subtracted dispersion 
relation. 

In order to calculate the discontinuous part of X™^(/c 2 ; m 1; m 2 , m 3 )n, we substitute a 
dispersion relation for il^- 2 \{k — p) 2 ; m 2 , m 3 ): 

iH^k- P r ;mi ,m 3 ) = i f d^ ^y '?■'"■» , (V.55) 

into Eq. ()V.53j) . Then, we obtain 

1 r°° 

IZ(k 2 ; mi ,m 2 ,m 3 ) u = ~ / dM 2 Im^ 2 )(M 2 ; m 2 . /,/, ) 

d 4 p 



' (m,2+m3) 2 

We take the discontinuous part of this equation: 



(27r) 4 p 2 — m 2 + it] (k — p) 2 — M 2 + irj 

dM 2 lmiI (2 \M 2 ; m 2 , m 3 )/ (2) (A; 2 ; mi, M). (V.56) 



oo 



1 r°° 

Imil^(k 2 ; m 1 ,m 2 ,m 3 ) u = - dM 2 lmiI i2) (M 2 ;m 2 ,m 3 )lmiI^ 2) (k 2 ;m u M). 

(V.57) 
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Using the expression for Irm/( 2 )(p 2 ; mm 2 ) we can easily evaluate Eq. (jV.57|) numerically. 

Next we turn to the continuous part. The continuous part of Eq. ()V.53j) is divergent. 
Therefore one needs corresponding counter terms in the Lagrangian. The second and third 
diagrams of Fig0] appear at two- loop order, which cancel the divergences in the bare sunset 
diagram, the first diagram of Fig0] In this paper, we do not explicitly go through the 
renormalization procedure but concentrate on the finite part. 

In order to calculate the finite part, we use the twice-subtracted dispersion relation: 



ReiT^(k 2 ; m 1 ,m 2 , m 3 ] 



11 



= Re <^ iTZHk ; m x ,m 2 , m 3 )n - «X S V ^(0; mi, m 2 , m 3 )u - k 



_d_ 

dk 1 



iT s ^(k ;mi,m 2 ,m 3 ] 
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k 2 =0 

7 (mi+m2+m3)2 M\M 2 - k 2 ) 1 • 

The second and the third terms in the left hand side are divergent and have to be renor- 
malized but the right hand side is finite and remain unchanged by renormalization. Using 
the results for Im£Z^£(M 2 ; mi, m 2 , m 3 ) we can easily evaluate Eq. (|V.58|) numerically. Af- 
ter Rei2^(m^ s ; mi, m 2 , m 3 )n is subtracted, Eq. (|V.58J) coincides with that in modified 
minimal subtraction scheme up to an irrelevant overall factor. 



B. T-dependent part 

We rewrite the integrand of the first term in the T-dependent part Eq. (|V.54|) so that it 
becomes the same form as that in the second term, i.e. (a delta function) x (some functions). 
This can be done by recombining the two factors in the integrand of F^ 2 '(k—p; m 2 , m 3 ) + (2 <-> 
3) and F^ 3 \k — p;m 2 ,m 3 ), T-independent part and T-dependent part of iAf x . The result 
is as follows: 

2^ n (fc;mi,m 2 ,m 3 )ii 

{n(p )27i5(p 2 -m 2 ) [l {2 \k - p-m 2 ,m 3 ) + (F {2 \k - p;m 2 ,m 3 ) + (2 <-> 3)) 



(2vr) 4 

+F (3) (fc -p;m 2 ,m 3 )] 

+n(po)27r5(p 2 — m 2 ,) [F 2 '(k — p; mi, m 3 ) + F^ 2 \k — p; mi, m 3 )] 

+n(p )2n5(p 2 - m 2 )/ (2) (k - p; m u m 2 ) } . (V.59) 
When we put k = (k , 0), Eq. (jV.59|) is reduced to 
T^ n (k , 0; mi, m 2 ,m 3 ) n 
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— y~] I / dun(u)^uj 2 - m\ I {2 \k Q + too, \Ju 2 - m?;m 2 ,m 3 ) 

^ T __j_ L^mi L 

+ (V 2 ) (/c + to, yju 2 -ml; m 2 , m 3 ) + (2 <-> 3)^ 
o + to, \]u 2 — m\; m 2 , m 3 ) 
do;n(o;) \J uj 2 — m\ I^ 2 \k + tuj, \Juj 2 — m\; mi,m 3 ) 



;r(2 
(2), 



+F ( } (A;o + to, a/cj 2 - m2;mi,m 3 ) 



+ 



cLm(u;) \J uj 2 — m| • I^ 2 \k + tuj, \J uj 2 — m§; mi, 7712) i , 



(V.60) 



Using the expressions of J (2) , F( 2 ) and F^ 3 ) obtained in the previous section, we can evaluate 
the continuous and discontinuous parts of Eq. (jV.60|) numerically. 

VI. CONTRIBUTION OF THE THERMAL SUNSET DIAGRAM TO THE SPEC- 
TRAL FUNCTION OF THE SIGMA MESON AT FINITE TEMPERATURE 

The purpose of this section is to see how two-loop diagrams affect observables by eval- 
uating the contribution of the thermal sunset diagram to o spectral function at finite tem- 
perature in the 0(4) linear sigma model. 

It is known that naive perturbation theory breaks down at T 7^ and that resummation of 
higher orders is necessary |23l24j . We adopt here a resummation technique called optimized 



perturbation theory (OPT) £Q. We first briefly review the procedure of OPT applied to 
0(4) linear sigma model. The original linear sigma model Lagrangian is as follows: 

(VI.61) 



^ = \ K 9 m&) 2 ~ A*V<] ~ T\^ 2 + h( ^° + counter terms, 



4! 

where fa = (a, tt) and hfa being the explicitly symmetry breaking term. For the renormalized 
couplings fi 2 , A and h and the renormalization point re we use the values determined in [Q: 
/J, 2 = -(283 MeV) 2 , A = 73.0, h = (123 MeV) 3 , re = 255 MeV. 

In OPT one adds and subtracts a new mass term with the mass m to the Lagrangian. 
Thus, we have 

£ = \{{d^faf - m 2 fal] + l -x4> 2 - ^{fa 2 ) 2 + hfa + (counter term), (VI.62) 

where x = m2 ~ ^ ■ The idea of OPT is reorganization of perturbation theory: one treats 
the added one as a tree-level mass term while the subtracted one as perturbation. 
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When the spontaneous symmetry breaking takes place, tree level masses of ir and a read, 
respectively, 

ml = m 2 + ^ 2 , ml = m 2 + ^ 2 , (VL63) 

where £ is the vacuum expectation value of o and determined by the stationary condition 
for the thermal effective action V(£,T, m 2 ) £Q: 

«^=a_a (VI.64) 

Note that the derivative with respect to £ does not act on m? . 

If Green's functions are calculated in all orders in OPT, they should not depend on the 
arbitrary mass, m. However, if one truncates perturbation series at a certain order they 
depend on it. One can determine this arbitrary parameter so that the correction terms are 
as small as possible. We adopt the following condition [T]: 

U^k 2 =ml) + U w (k = 0;T) = 0, (VI.65) 

where the first and second terms are respectively T-independent part and T-dependent part 
of the one- loop self-energy of 7r. 

Let us now turn to the discussion on the spectral function of a defined by Eq. ()H.24j) : 

(h V\= -- ImII CT (A;o,k) 

Here n CT (fco, k) is the self-energy of a. Its real and imaginary parts are related with (1,1) 
component, n^. 1 , via Eqs. flII.13|) and (|II.14|) . 

As was already mentioned, the spectral function at one-loop order was studied by Chiku 
and Hatsuda £Q . We want to see how the spectral function at one- loop order is modified by 
adding the thermal sunset diagram. Thus, we take 

n^o, k) = n^fo, k) x _i oop + Ul\kQ, k) sun . (VI.67) 

The first term is the renormalized one-loop self-energy calculated in P . The second term is 
the renormalized sunset diagram depicted in Fig|S]and given by 

A 2 

n"(/c ,k) sun = -—iT sun (ko,\i; 171^,1710-,:, m 0a )ii- (VI.68) 
We show p a (k ,k = 0) at T = 200 MeV and T = 145 MeV in Fig© At T = 200 MeV, 
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the spectral function at one-loop order consists of a 5-function peak for a and a continuum. 
By including the sunset diagram a acquires a width of the order of 10 MeV. At lower 
temperature, T = 145 MeV, an enhancement of the spectrum near the threshold is observed 
at one- loop order. When we include the sunset diagram, this feature is not lost. 

Let us discuss the above results. At high temperature (T = 200 MeV), the mass of a is 
smaller than twice the pion mass and the decay a — > tttc is forbidden. As a result, within the 
one-loop calculation o has zero width. However, at finite temperature a can interact with 
thermal particles in heat bath and change into other states. Among such processes, those 
which are taken into account by including the sunset diagram are represented by Fig|B]with 
($,01,02,03) assigned to, for example, (a, a, 7r, tt). The processes which correspond to (a) 
and (b) in Fig|3]are possible for k > 2m 07r + ni 0a . This affect the spectrum at high energy, 
(c) and (d) in Fig|3]drop off for positive k . (e) and (f) affect the spectrum at low energy 
since they are possible for fc < m 0a ~ ^ m oir- The processes which correspond to (g) and (h) 
in FiglHl are shown in Fig0 They are the most important and give a finite width to a since 
they are allowed at arbitrary positive fco- However, we observe that their effects at lower 
temperature (T = 145 MeV) are small. The reason is traced back to Eq. (jIV.50j) . The term 
representing the probabilities for (g) and (h) in Fig|3] is the fourth term. That integral at 
lower temperature is suppressed due to the statistical weight since at lower temperature the 
masses of a and 7r are large. 

Finally we note that, as a consequence of the non- vanishing discontinuous part of the 
sunset diagram, the spectral function is also non- vanishing in the all range of k$. 

VII. SUMMARY AND PERSPECTIVE 

We have studied the sunset diagram for scalar field theories at finite temperature in the 
real time formalism. We have explained how we can reduce it to an expression written 
in terms of one-loop self-energy integrals, which can be easily evaluated numerically. We 
have also discussed what physical processes are contained in the discontinuous part of the 
diagram. We have found that there exist processes which occur only at finite temperature 
and some of them are allowed at arbitrary energy. As a result, the discontinuous part of the 
sunset diagram is non-vanishing in all the energy region, which is a remarkable feature at 
finite temperature and manifests itself at two-loop order. 
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As an application of the result, we have demonstrated how the spectral function of a at 
finite T at one-loop order is modified when we include the thermal sunset diagram. At high 
temperature, where o — > ttit is forbidden, a acquires a finite width of the order of 10 MeV 
due to collisions with thermal particles in the heat bath while a does not have a width at 
one-loop order. At lower temperature the spectrum is almost unchanged. 

Finally we comment on the effect of other two-loop diagrams on the spectral function. We 
have seen that the threshold enhancement, which was first found in the one-loop calculation, 
is retained if we include the sunset diagram. However, in the present calculation the effect 
of the thermal width of 7r in the a — > Tin decay is not included. This is taken into account 
by including the diagrams such as shown in FiglHl i n which an internal 7r changes into o 
absorbing a thermal 7r in the heat bath. 

In the one-loop calculation, tt has a width, = 50 MeV, at the temperature, T = 
145 MeV, at which the threshold enhancement is observed for a. If we include this effect as 
a constant complex mass shift for 7r in the one-loop self-energy for a, we expect a to acquire 
a width twice as that for it, i.e. r CT = 2T n |25j . This implies that when the two- loop diagrams 
such as FiglHlare included the spectral function for a would be significantly modified with 
the width of about 100 MeV. The calculations of those diagrams are now in progress [2*6*] . 
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FIG. 3: The amplitudes in Eq. (|IV.50|) responsible for the disappearance and reappearance of 
(f>i stands for an anti-particle of (pi. (a) minus (b) corresponds to the first term in Eg. (|lV,50|) . (c) 
minus (d) to the second, (e) minus (f) to the third and (g) minus (h) to the forth. Amplitudes 
obtained by permuting particle labels of (e), (f), (g) and (h) also exist. 
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FIG. 4: Renormalized T-independent part of the sunset diagram. The second and the third 
diagrams cancel the sub- and overall divergences in the first diagram (bare sunset diagram) respec- 
tively. 




FIG. 5: The sunset diagram for a. Solid and dashed lines correspond to a and n respectively. 
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FIG. 6: Spectral function of a p a {ko,k = 0) at T = 200 MeV (upper panel) and T = 145 MeV 
(lower panel). Solid line corresponds to p a at one-loop order and dashed line to that with one-loop 
self-energy and the sunset diagram. 
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(a) (b) (c) 

FIG. 7: The processes which are allowed at all positive ko contained in the discontinuous part of 
the sunset diagram for a (Fig|5J). 




FIG. 8: two-loop self-energy diagram in which internal ir changes into a by absorbing thermal tt 
in the heat bath. 
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APPENDIX A: CALCULATION OF Im^ 2 )(p;mi,m 2 ) AND F^\p;m 1 ,m 2 ) 



In this appendix, we derive Eas. (llll.35j) - ()111.37j) and Eas. fjlll.39D - fllll.41D . 



From Eq. (jlll.28|) we obtain 



1 f d k 

lmiF {2) (p;m u m 2 ) = - / — — ^^({p + k) 2 - m\)n(k Q )2Ti5(k 2 - m\) 




x / dcos95(p 2 + 2p A;o-2|p||k|cos(9 + k 2 -ml). (A.l) 



/ dk d\k\\k\ 2 n(k )5(k 2 - m 2 2 ) 



ImiF^ (p; mi, m 2 ) receives contribution from k which satisfies 

\p 2 + 2p k + m 2 , - ml\ 

2isiki < 

and k 2 , = k 2 + m\. We make a square of Eg. f)A.2j) and obtain 



(A.2) 



4po^o(p 2 + ml — Tnf) < —{p 2 + ml — m\) 2 — 4p 2 1 m 2 ! — 4p 2 |k 



(A.3) 



There are three cases when Eg. (jA.3J) holds: 



A. LHS> and RHS> 0, (RHS) 2 -(LHS) 2 > 0. 

B. LHS< and RHS< 0, (RHS) 2 -(LHS) 2 < 0. 

C. LHS< and RHS> 0. 

First, we calculate (RHS) 2 -(LHS) 2 . 
(RHS) 2 - (LHS) 2 

= 16p 4 |k| 4 + {8p 2 [(p 2 + m 2 2 - ml) 2 + 4p 2 m 2 ] - 16(p 2 + m 2 - m 2 ) 2 p 2 ] |k| 2 

, 17 2 , 2 2\2 , a 2 21 2 irf 2 , 2 2\2 2 2 

+ [(p + m 2 — mjj + Ap m 2 \ — lb(p + m 2 —m l )p Q m 2 
= 16p 4 |k| 4 + 8 [(p 2 - 2p 2 )(p 2 + m 2 2 - ml) 2 + Ap 2 p 2 m 2 2 ] \k\ 2 

+ [(p 2 + m 2 2 - ml) 2 - Ap 2 m 2 2 ] 2 (A.4) 



Whether there exists |k| 2 which satisfies (RHS) 2 — (LHS) 2 < depends on the sign of the 
following expression: 



D — 16 [Qo 2 - 2p 2 ){p 2 + m 2 2 - ml) 2 + Ap 2 p 2 m 2 2 ] 2 - 16p 4 [{p 2 + m 2 - m 2 ) 2 - 4p 2 m 2 ,] 2 



64po|p| 2 G° 2 + m l - m l) 2 [p 2 - ( m i + m 2) 2 ] [p 2 - (mi - m 2 ) 2 ] . 



(A.5) 
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Thus, for (mi — m 2 ) 2 < p 2 < (mi + m 2 ) ; 



(RHS) 2 - (LHS) 2 > 0, 



for p 2 < (mi — m 2 ) 2 or (mi + m 2 ) 2 < p 2 



|k| 2 <|k| 2 or |k| 2 > (RHS) 2 - (LHS) 2 > 0, 

|k| 2 < |k| 2 < \k\ 2 + <^> (RHS) 2 - (LHS) 2 < 0, 



where |k| + are given by 



|k|l 




1 - 



(m 2 + m x ) 2 

p2 



1 - 



(mi — m 2 )' 
p2 



\l 



Secondly, for p 2 > m\ — m 2 . 



p k Q ^ 



LHS ^ 0, 



for p 2 < m 2 — m\ 



Thirdly, for p 2 > 



LHS ^ 0. 



RHS < 



for p 2 < 



| k |2 ^ (p 2 + m 2 - m 2 ) 2 + 4p 2 m 2 = |k|2 



RHS ^ 



-Ap 2 ~ |H0 

Therefore, the conditions for the above three cases are respectively given by 

A. p 2 < 0, Pok (p 2 + m 2 , — m\) > and |k| 2 > |k|+. 

B. p 2 > (mi + m 2 ) 2 , p ko(p 2 + m 2 , - m\) < and |k| 2 < |k| 2 < |k| 2 , 

< p 2 < (mi — m 2 ) 2 and Pok (p 2 + m\ — m\) < and |k| 2 < |k| 2 < |k| 2 , 
p 2 < 0, p k (p 2 + m 2 , - m\) < and |k| 2 < |k| 2 < \k\ 2 . 

C. p 2 < 0, p k (p 2 + m 2 - m 2 ) < and |k| 2 > |k|„. 

The case C can be combined with the third of the case B as 
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p 2 < 0, p k (p 2 + m 2 . - ml) < and |k| 2 > \k\ 2 _. 

1. For p 2 > (mi + m 2 ) 2 or < p 2 < (mi — m 2 ) 2 , 

1 /" |k| + 

IrmF' 2 ** (p; mi, m 2 ) 



8t|P| J\k\- 
1 

dwn(cj) 



87r|p| 

1 1 

16tt|p|^ 



In 



1 - e -/3t 



1 - e-^- 



2. For (mi — m 2 ) < p < (m x + m 2 ) , 



lmiF^ 2 \p;m 1 ,m 2 ) = 0. 



3. For p 2 < 0, 



ImiF^ (p; mi, m 2 ) 



1 



8tt| P | 
1 

Nil 
i 



|k|- 

OO /"OO 



rf|k| — n(u k ) 
|k|+/ ^fc 



dcun(c<j) 



16tt| P | /3 



■ln|(l-e-^+)(l-e-^-)| 



where 



Similarly, 



mj — m? \ 2 
1 + bo|± 



1 - 



(m 2 + mi) 5 

^2 



1 - 



(mi - m 2 y 
p 2 



1. For p 2 > (m x + m 2 ) 2 , 

IrmF® (p; m 1 ,m 2 ) 



1 /"' k ' + Ikl 

8vr|p| J| k |_ ^fc 
1 

dcjn(cj)n(|p | — 



8tt| P 
1 



1 



87r|p|/3e0l*>l - 1 
2. For (mi — m 2 ) 2 < p 2 < (m x + m 2 ) 2 , 



In 



X _ e -/3^+ e /?(|po|-w_) _ x 



1 — e -/3^- e /3(|po|-w+) _ I 



lmiF^\p;m 1 ,m 2 ) = 0. 
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3. For < p 2 < (mi — "m 2 ) 2 , 
ImiF^ (p; mi, 777.2) 



8tt| P | 
1 

8tt| P | 
1 



|k|+ | k | 

d\k\ — n(uj k )n(\p \ + u k ) 

|k|_ 



U)+ 



dum(uj)n(\po\ + cj) 
1 



In 



/3 - 


l 


1 - e~^+ 


- e^ln 


1 - e~^- 



X _ e -/9(bo|+w+) 



1 _ e -/9(|p |+w_) 



. (A.13) 



4. For p 2 < 0, 



ImiF^lp; mi, ms) 

1 / Z" 00 Ikl Z" 00 Ikl 

— — - / d|k| — n(u k )n(\p \ + u k ) + / d|k| — n(o; fc )n(|p | - u k ) 
8tt|p| VJ|k|_ ^fc i|k|+ ^fc 

/>oo 

dum(u;)n(|po| + uj) + / dum(u;)n([po| — uj) 

J w+ 

1 



8tt| P | 
1 



1 

4 



+ e -/3bol ln | 1 _ e - 



/3(w+-|po|)| 



o/3|pol 



— [-ln|l 



+ e /3|pol ln|l -e 



-/3(w_+|po|) 



I] ■ 



(A. 14) 



APPENDIX B: PHYSICAL MEANING OF |k|± 



In this appendix we explain the physical meaning of |k|± defined by Eq. (jA.6l) . 
Let us first consider the case p 2 > (mi + m<i) 2 or < p 2 < (mi — mj) 2 . For simplicity we 
suppose po > 0. Since internal particles are on shell, in the center-of-mass frame we have 

2 



P 



m; 



|k| 2 ± \/ml 



where + is for p 2 > (mi + m2) 2 and — is for < p 2 < (jn% — m 2 ) 2 . In either case 

2 _ [p 2 - (mi + m 2 ) 2 ][p 2 - (mi - m 2 ) 2 ] 
1 1 ~ 4^2 • 



(B.l) 



(B.2) 



If we boost the system to the positive z-direction, the external and internal momenta become 



p' = p cosh6>, 

P± = 0, p' z = p smh.6 (9 > 0), 

k' = kocosh.6 + k z sinh6, 

k'j_ = kj_, k' — k z cosh8 + /cosinh^. 



(B.3) 
(B-4) 
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Since — Ikl < k 7 < Ikl 



where 



|k'| 2 < |kf < \k'\ 2 + , 



(B.5) 



|k'|| 



±|k|cosh# + |A; |sinh#)" 
±| k |^ + N ^ 2 



Po 



Po 



(mi + m 2 ) 
p2 



1 



(mi - m 2 ) 
P 2 



Po + 



1 + 



2 2 

m,2 — vti\ 



p 



IP 



(B.6) 



This is nothing but Eq. (|A.6|) . Therefore, for p 2 > (mi + m 2 ) 2 or < p 2 < (mi — m 2 ) 2 , |k| + 
and |k|_ are respectively the maximum and minimum values of |k| such that the internal 
particles are on-shell. 

Let us next consider the case p 2 < 0. In the Breit frame, p = (0, p), with p in the positive 
z-direction we write the momentum of the particle with mass m 2 by k = (k , k) where k can 
be positive or negative. Then, the momentum of the particle with mass p + k = (ko, p + k). 
Since internal particles are on shell, we have 



k 2 = m 2 2 + k 2 z + \k ± \ 2 = m\ + (p z + k z ) 2 + |kj 



(B.7) 



From this equation we obtain 

k z = 



p 2 + m\ — 777-2 — P 2 + m i ~ m 2 



2p z 



k 2 



-p 2 + (mi + m 2 ) 2 ][-p 2 + (mi - m 2 ) 2 ] 2 2 

=v +l ±l ' 



When boosting the frame to the z-direction, we suppose the momenta become 



(B.8) 



P = (Po, P')> k ' = ( fc 0) k ')> P +k' = (p + k' , p + k') 



(B.9) 



Then, p' and k' are given by 



p' Q = Iplsinh^, 
|p'| = Ipjcosh^, 

k' = fcocosh^ + fc^sinh^, 

k'j_ = kj_, k' z = k z cosh9 + kosinh.9. 



(B.10) 
(B.ll) 
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Using Eqs.(|R8}, (lRlO|) and (lB~TT| . we obtain 



k'| 2 = (k z cosh8 + fc sinh6>) + |kj 



(B.12) 



|k'| 2 has a minimum when | kj_ | 2 = but not a maximum. Depending on the sign of ko, the 
minimum value is given by 

2 



Ml 



k-. 



IP 



P 



±Jk 2 + m 2 - Po 



-p 



1 

4 



2 2 

m\ — m 2 



—p z 



|P|± 



1 + 



(mi + m 2 ) 

_p2 



1 + 



(mi — m 2 ) 2 
— p 2 



Po} -(B.13) 



Eq. (jB.13|) is nothing but Eq. ()A.6|) . Therefore, in the case of p 2 < 0, both of |k| + and |k| 
are the minimum values of |k| such that the internal particles can be on-shell. 
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